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Abstract 

The analysis of flowing afterglow plasmas containing H 3 + or D 3 + ions indicates 
that the de-ionization of such plasmas does not occur by simple dissociative 
recombination of ions with electrons. An alternative model of de-ionization is proposed 
in which electrons are captured into H 3 ** autoionizing Rydberg states that are stabilized 
by collisional mixing of the Rydberg molecules' angular momenta. The proposed 
mechanism would enable de-ionization to occur without the need for dissociative 
recombination by the mechanism of potential-surface crossings. 

I. Introduction 

The dissociative recombination (DR) of electrons with H 3 + ions 

H 3 + + e- H+H+H 

H 2 +H (1) 

is thought to be a key process in the ionospheres of the outer planets, in interstellar 
plasmas, and in many electrical discharges and its recombination coefficient enters as an 
important parameter in models of such plasmas. Since H 3 + is the simplest polyatomic 
ion, it has also become a crucial test case of the theory of DR. While considerable 



efforts have been devoted to this problem by theorists and experimentalists, a clear 
picture has not yet emerged and the experimental results still present a somewhat 
confusing picture. It is not even clear if binary recombination actually takes place. It was 
the basic goal of this work to examine the results of earlier experiments and to resolve 
some of the apparent paradoxes. In the course of doing this, it became necessary to 
repeat some earlier flowing-afterglow measurements. 

We will begin with a brief review of the history of experimental work on H 3 + 
recombination and its interplay with theory. A recent review by Dalgamo [1] outlines 
the role of H 3 + in terrestrial and extraterrestrial environments and his paper contains 
extensive references to earlier theoretical and experimental work on this ion. The brief 
"thumbnail sketch of H 3 + " by Oka and Jagod [2] provides an excellent introduction to the 
structure and vibrational states of H 3 + . A fairly complete account of the theory and vari- 
ous measurements of dissociative recombination rates can be found in a paper by Bates 
[3 ]. We have compiled a list of measured rate coefficients and related information in 
Table 1. 

From the point of view of theory, the problem is that the direct mechanism [4] of 
DR fails in the case of H 3 + ions. Calculations of the relevant potential surfaces [5] 
showed that there are no crossings between the ground-state ionic potential surface and 
repulsive states of the neutral molecule that would facilitate DR of H 3 + in v=0. The 
indirect mechanism, which is thought to involve an intermediate, autoionizing and 
predissociating Rydberg state, relies on the same crossings of potential surfaces and is 
subject to the same problem. Such arguments led Michels and Hobbs [5] to propose that 
H 3 + ions in v=0 should not recombine with electrons efficiently, but that vibrationally 
excited H 3 + ions in v>3 might do so. Similar predictions had been made earlier by 
Kulander and Guest [6]. In contrast, the only experimental data (obtained by using the 
microwave afterglow method) that were available at that time (Leu et al. [7] and 
Macdonald et al. [8]) gave recombination coefficients on the order of 2x10*7 cm-Vs, 



which is a fairly typical value for a small species. Since it appeared reasonable to assume 
that the experimental data referred to vibrationally unexcited H 3 + ions, the results were 
clearly in conflict with Michels and Hobbs' expectation. An apparent resolution of the 
dilemma was proposed by Adams et al.[9]. Their observations in flowing afterglow 
plasmas suggested that two distinguishable forms of H 3 + ions might exist, a major species 
that recombined fairly fast and a minor species that recombined only very slowly (a<2 
xlO'8 cm3/s). This observation permitted the interpretation that the slowly recombining 
ion fraction consisted of H 3 + (v=0) ions and it seemed natural then to identify the fast 
recombing fraction as H 3 + in vibrationally excited states in v>3. This interpretation 
gained acceptance in spite of the fact that Adams et al.[9] had provided no independent 
experimental evidence for the presence of vibrational states or for their identification. 

Several attempts were made to explain the discrepancy between the experimental 
results of the microwave-afterglow and the flowing-afterglow method. Since "impurity" 
ions are invariably present in afterglow plasmas, experimentalists examined their effects. 
Two problems can arise: Conversion of H 3 + to more complex ions tends to increase the 
rate of electron loss, while the presence of non-recombining ions (e.g. He + , H + ) would 
reduce the electron loss. Since ions of mass 17 amu (most likely CH 5 + ions) were known 
to be present in the microwave afterglow experiment of Macdonald et al.[8], Johnsen [10] 
mentioned the possibility that those measurements could have been affected by impurities 
without giving an estimate of the magnitude of the effect. A reanalysis of those data 
which includes the effect of CH 5 + recombination, indicates that Macdonald et al.'s value 
of a(240K) should be corrected downwards by only about 10%. The impurity effect 
would have been larger at the higher electron temperatures. 

The second "impurity" problem has been discussed by Canosa et al.[ 1 1,12] and 
by Smith and Spanel[ 13]. It is one peculiarity of He/H 2 afterglow plasmas that He + ions 
are not efficiently removed by ion-molecule reactions, since the reaction of He'*' with H 2 
is unusually slow (k= 1.1 x lO' 1 ^ cm3/s)[14]. Therefore, non-recombining He + ions can 
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accumulate in the late afterglow and reduce the recombination rate. It appears that the 
original experiments by Adams et al. [ 9] were somewhat affected by this problem, but 
later work by Smith and Spanel [13] showed that this is not the only reason why the 
recombination rate is reduced in the late afterglow. 

The conjecture was put forth by Johnsen [unpublished, but a brief reference sur- 
vives in the paper by Canosa et al.[l 1] that perhaps the recombination of H3 + with elec- 
trons might produce H3* molecules in highly excited Rydberg states and that these 
molecules would be field-ionized by the Langmuir probes that are used in flowing 
afterglow experiment to measure electron densities. As a consequence, the current 
collected by the probes might be erroneously large and the recombination loss would be 
appear to be reduced. The conjecture, however, cannot quantitatively account for the 
observed effects since the number of electrons that can be collected from Rydberg 
molecules is limited by the diffusive particle current to the probe surface. Due to the 
larger mass of H3* molecules, their thermal flux is far smaller than that of free electrons 
(by the factor [m e /m(H3 + )] , ^ = 0.04). The conjecture, nevertheless, may have contained a 
kernel of truth. We will return to this subject in a later section (V.2). 

The first clear evidence that Adams' et al. [9] identification of H3 + vibrational 
states could not be upheld came when Amano [15] published results of an infrared 
absorption study in a decaying hydrogen plasma. He concluded that H3 + ions in the 
spectroscopically- identified vibrational ground state recombined with electrons about as 
fast as had been measured in the microwave experiments. Several objections to Amano's 
results were raised, for instance that the electron density in his plasmas was sufficiently 
high to enable collisional-radiative recombination. A further objection was that Amano's 
assumption of equal electron and H3 + densities might not have been justified, since the 
plasma could have contained ions other then H3 + . Amano addressed the questions in a 
second paper 16 and pointed out that the ion density followed a hyperbolic decay (electron 
density proportional to reciprocal time) very well over a large range. This would not be 



expected, if either of the two objections was valid. The absence of collisional radiative 
recombination, however, presents a puzzle. Bates et al. [3] have discussed this question in 
more detail. A further infrared absorption experiment in hydrogen plasmas has recently 
been carried out by Feher et al. [17] . Those authors only concluded that Amano's 
recombination rate should be considered as a upper limit which is not much of a 
constraint. 

Additional measurements of the recombination of H 3 + in flowing-afterglow 
plasmas were carried out by Canosa et al.[l 1,12 ] Their technique was rather similar to 
that of Adams et al. [9] but it differed in the method of data analysis. The measured 
recombination coefficient were close to those obtained by Amano [15,16] and by the 
early microwave afterglow experiments. 

Several measurements of recombination cross sections for H 3 + were made by 
using the merged-beam technique [18, 19, 20, 21 ]. The outstanding virtue of this 
technique is that the energy resolution (in the electron-ion center-of-mass system) is 
higher than that of plasma afterglow experiments and that fine structures in cross sections 
become amenable to study. Unfortunately, the only control of the internal energy of the 
reacting ions (e.g. vibrational energies) consists of modifying conditions in the ion 
source. The results of Hus et al.[l 8] showed that the recombination cross section of H 3 + 
was strongly affected by the conditions in the ion source and this was ascribed to 
vibrational excitation. In some measurements made with the merged-beam method 
(Yousif et al.[21]), the vibrational state of the recombining ions was deduced from the 
threshold energy for electron-ion dissociative excitation. The result seemed to indicate 
that the H 3 + ions were in v= 0 and that those recombined only slowly with electrons. 

A major new innovation in experimental techniques occurred when ion storage 
rings [22, 23] became available for measurements of dissociative recombination. In these 
large-scale instruments, the residence time of ions can be made sufficiently long (8 
seconds) to allow vibrationally excited ions to relax to the vibrational ground state by 
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radiative transitions. The authors do not explicitly state this, but presumably even the 
vj=l metastable level with a calculated lifetime (Dinelli et al.[24]) of 1.2 s would have 
decayed by infrared transition to the V2=l level. The work by Larsson et al. [22] and by 
Sundstrom et al.[23] using the CRYRING storage ring yielded recombination cross 
sections for H 3 + (v=0) over the wide energy range from 0.0025 to 30 eV. The inferred 
recombination coefficient of (1.15 ± 0.13) xlO'^ cm-Vs at 300 K (as calculated by 
Sundstrom et al. [23]) again falls in the range of the higher values. At present, the energy 
resolution in this technique is limited by the electron motions perpendicular to the beam. 
The effective "temperatures" (kT = 0.1-0.15 eV) are actually higher than in the afterglow 
experiments and it should be kept in mind that the conversion from measured cross 
sections to thermal rate coefficients requires an extrapolation of the cross sections to low 
energies. 

At this time, it appeared that the problem had been resolved in favor of the larger 
recombination rates. The conflict between theory and experiment, sometimes referred to 
as an "enigma" [3], appeared to fade into the background when new recombination 
mechanisms [25, 26] were proposed that seemed capable of accommodating higher 
recombination rates for H 3 + (v=0 ) ions. The situation changed again when Smith and 
Spanel [13] presented new data in support of Adams et al.s' [9] contention that only H 3 + 
ions in v>3 recombined efficiently with electrons and that v=0 would recombine only 
slowly. Their experiments included many corroborating tests and the conclusions seemed 
to be well-founded. Our work was, in part, motivated by that experiment. 

In the article, we will present new results of flow-tube measurements on the 
recombination of H 3 + and D 3 + ions. While we were able to reproduce some of Smith 
and Spanel's [13] experimental results, we will show that their interpretation of the 
anomalous decay of H 3 + plasmas cannot be correct. We were eventually led to 
conclude that the peculiar decay of H 3 + plasmas probably has little to do with vibrational 
excitation of the ions. A different and more interesting explanation will be offered. 



In the following, we will make a distinction between the customary 
"recombination coefficient" of a specified ion (denoted by "a") and the "de-ionization 
coefficient" (denoted by "p”) of a plasma. The de-ionization coefficient is meant to 
represent the total measured volume loss (excluding diffusion) of free electrons. It may 
refer to a mixture of ions and it may include other processes, for instance re-ionization of 
recombination products. This distinction does not bias the discussion in favor of a 
postulated recombination process, and it avoids the use of ill-defined concepts, such as a 
time-dependent recombination coefficients. It will be seen later that this distinction is 
important in the case of H 3 + plasmas. 

II. Experimental Methods 

The flowing afterglow method to determine electron-ion recombination coeffi- 
cients and products of recombination is quite well developed [9], The emphasis in this 
paper is on the anomalous features of the decay of H 3 + plasmas. Therefore, we will 
devote more space to questions of vibrational excitation of H 3 + ions and the methods of 
data analysis. 

Our flow tube (see Fig. 1) is a stainless steel tube (diameter 6 cm, length = 36 
cm ) with two attached long quartz windows for spectroscopic observations (not used in 
the present experiment). As is customary, reagent gases are added through several 
movable and fixed reagent inlets. A movable Langmuir probe (diameter =25|im, length = 

0.28 cm) serves to measure electron densities. We have given a detailed discussion of the 
use of Langmuir probes in flow tubes in a recent publication [27]. A quadrupole ion 
mass spectrometer samples plasma ions from the downstream end of the flow tube. All 
gas pressures, flow rates were measured by capacitance manometers and calibrated gas 
flow meters. The carrier gas was helium (at pressures near 1 Torr). Flow velocities 
ranged from 5000 cm/s to 6000 cm/s. The gas temperature was typically 295 K. 



II. 1 Method of plasma generation and ion conversion 

A plasma in the pure carrier gas helium is generated in an upstream microwave 
discharge. Argon gas (at a pressure of typically 15 % of the helium pressure) is added at 
a downstream gas inlet to convert metastable helium atoms to argon ions by Penning 
ionization and to reduce the rate of ambipolar diffusion. In the region downstream from 
the argon inlet, the plasma then contains mainly thermalized electrons, Ar 1- ions and some 
He + and He 2 + ions (< 0. 1 % of Ar 1- , as measured by the mass spectrometer). Such a 
plasma is essentially non-recombining. At the position of the reagent inlet, the electron 
density was typically 4xl0 10 cm' 3 . 

To create H 3 + ions, H 2 is added at the reagent inlet at densities of 10 14 cm' 3 . 
The following reaction sequence occurs: 

Ar + + H 2 — > Ar H + + H (+ 1.53eV) (2) 

ArH + +H 2 -> Ar + H 3 + (+ 0.57 eV). (3) 

The exoergicities refer to Ar 4- in the 2 P 3 /2 fine structure state and ground vibrational and 
electronic states for the molecular species. Thermochemical data were taken from Radzig 
and Smirnov [28] . The rate coefficients for reactions 2 and 3 are on the order of 10'^ 
cm 3 /s [29], 

II.2 Vibrational excitation and quenching 

The question of vibrational excitation in the H 3 + ions needs to be examined in 
detail. If it is assumed, which however seems unlikely, that the entire energy of 2.1 eV 
that is released in the two-step process appears in vibrational excitation of the H 3 + 
product, then H 3 + could be produced in vibrational states up to v=5. However, the more 
highly excited H 3 + ions will be destroyed quickly by proton transfer to Ar, 
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H 3 + + Ar -> Ar H+ + H 2 (-0.57 eV), (4) 

which becomes energetically possible for H 3 + ions with internal energies above 0.57 eV. 
There are only three H 3 + vibrational states below 0.57 eV, the ground state [Ai (0, 0°)], 
the v 2 =l bending-mode vibration [E (0, l 1 ) at 0.3126 eV], and the vi=l breathing mode 
vibration [Ai (1,0°) at 0.394 eV] (see Lie and Frye [30 ] or Oka [2]). The corresponding 
energies for D 3 + are 0.2275 eV (v 2 =l) and 0.2856 eV (vi=l). 

Since the density of argon in the experiment is typically 5x 10 15 cnr^, and exo- 
ergic proton transfer is invariably fast (k ~10' 9 cm 3 /s), destruction of the higher vibra- 
tional states would occur in a time of less than 1 psec. This time is far shorter than the 
time over which recombination occurs. 

One may wonder, however, if vibrationally excited H 3 + ( v>l ) might enter a reac- 
tion equilibrium with vibrationally excited ArH + ions (via reaction 3 in both directions) 
and thus would be constantly regenerated. It seems very unlikely that vibrational excita- 
tion would persist very long in such a process since in each proton transfer some vibra- 
tional energy will be converted to translational energy. Furthermore, ArH + ions should 
then be detected by the mass spectrometer, but there is no experimental evidence for a 
significant presence of ArH + ions. The conclusion is that H 3 + ions in v>l can excluded 
from further consideration. 

The radiative lifetime of the vi=l level is very long since a direct transition to the 
ground state is dipole forbidden. Dinelli et al. [24] calculated a lifetime of 1.2 s for 
cascading to the V 2 = 1 level which then decays in 4 ms to the ground state. Thus, radiative 
decay of the vi=l level will not occur on the time scale of recombination (it takes about 1 
ms for the electron density to decay from 4xl0 10 cm* 3 to 1/8 of that value if the 
recombination coefficient is ot=2x!0' 7 cm^/s). Radiative decay of ions in the v 2 =l level 



is faster (x = 4 ms), but unless quenching occurs some could survive into the late 
afterglow. 

Collisional destruction (quenching) of v=l ions in collisions with H 2 will be 
considered next. Amano's [31] absorption studies of the V 2 hot band indicate that the 
V 2 =l level is quenched rapidly, with a rate coefficient of approximately 3xlO* 10 cm 3 /s. 
This eliminates ions in this state from further consideration, since they will be destroyed 
in 30 [is at [H 2 ]= lxlO 14 cm* 3 . On the other hand, there is considerable evidence (see 
Lee et al. [32], Blakey et al. [33] , and Smith and Spanel [13]) that some vibrationally 
excited ions are quenched with a far smaller rate coefficient of roughly 10*12 cm 3 /s. 
Those ions would survive for 10 ms at [H 2 ] = 1 xlO 14 (1 ms at [H 2 ]=l xlO 15 cm* 3 ) 
and could be present during the recombination phase of the afterglow. It seems very 
likely that the slowly-quenched ions are those in the vi=l level. A similar conclusion 
was reached by Bawendi et al. [34 ] who studied the absorption of H 3 + hot bands in 
He/H 2 discharges. They ascribed the slow quenching of the vj=l state to the lack of a 
dipole moment for transition to the ground state. 

The foregoing discussion shows that H 3 + (vi=l) ions may survive long enough to 
undergo recombination. An estimate of their relative abundance clearly would be of in- 
terest. The for the present purpose most relevant estimate is that of Lee et al.[32] . It 
refers to a flowing afterglow plasma of the same gas composition as that in our experi- 
ment, but of lower electron density. Their studies of ion-molecule reactions of H 3 + 
indicated that 70% of the ions were in the vibrational ground state and 30% in an 
(unidentified by them) excited state. Since the V 2 =l ions are quenched rapidly, we 
suggest that the 30% long-lived ions were in the vi=l state. 

A further, but seldom mentioned vibrational quenching may occur by superelastic 
collisions with electrons 


H 3 + (v=1) + e* — » H 3 + (v=0) - 1 - e* + kinetic energy . 


(5) 



The reaction would involve the same intermediate H 3 ** compound states that are 
invoked in indirect dissociative recombination, but it would be terminated by 
autoionization, rather than by dissociation. Unfortunately, nothing appears to be known 
for this type of quenching in the case of H 3 + , but Guberman [35] has carried out 
calculations for the corresponding process of N 2 + (v=l) ions. The surprising result was 
that the quenching rate exceeded the rate of dissociative recombination by a factor of 18. 
If the same were true for H 3 + (v=l) then most v=l ions would be quenched before 
suffering recombination. It is interesting to note that v=l ions account for only about 4 to 
6 % of all ions in He/H 2 discharges [34], Thus, the abundance of 30% given by Lee et 
al.[32] may be an upper limit that refers to plasmas of low electron densities where 
electron quenching is unimportant. 

In their discussion of the same reaction sequence. Smith and Spanel [13] did not 
consider the destruction of the higher vibrational levels via reaction (4). They concluded 
that most H 3 + ions would be highly vibrationally (v>3) excited. We reach the opposite 
conclusion that only a few (at most 30%) might be excited to vi=l. Experimental tests 
pertaining to vibrational excitation will also be discussed in Sec. III. 

II. 3 Para and ortho H 3 + . 

There is no a priori reason to expect that the para (three proton spins aligned) and 
the ortho (two proton spins aligned) forms of H 3 + should have the same recombination 
coefficients. In an afterglow experiment, two distinct recombination coefficients would 
be observable only if the para/ortho equilibration time is longer than the time scale of re- 
combination. In the opposite case of rapid equilibration, the concentration-weighted av- 
erage of the two coefficients would be measured. According to Oka [36], the conversion 
from one form into the other by proton exchange with H 2 , 
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H 3 + (ortho) + H 2 o H 2 + H 3 + (para) (6) 

is a rapid process, but unfortunately no rate coefficient is given. We assume that the pro- 
ton-transfer rate is sufficiently large (lolO' 1 ^ cm 3 /s would be sufficient at [H 2 ] = 10 14 
cm -3 ) to maintain the para-to-ortho ratio at its 1:1 equilibrium value. 

II.4 Presence of Hs + ions 

There is no reason to be concerned about formation of H 5 + ions at room 
temperature since the equilibrium constant for three-body association 

H 3 + + H 2 + He H 5 + + He (7) 

has the small value of 2xl0' 19 cm 3 at 300 K [37]. At the highest H 2 concentrations used 
in our experiments (1.7 xlO^ cm -3 ) the H 5 + /H 3 + ratio would be only 3xl0‘ 4 . 
Vibrationally excited H 3 + ions should not associate with H 2 , even at reduced 
temperatures, since the binding energy of H 5 + ( 0.35 eV) [37] is very close to the 
vibrational quantum in H 3 + . 

II. 5 Impurity estimates 

Gaseous impurities can be problem in afterglow experiments. Their effect is 
more severe in studies of H 3 + ions since the ions readily transfer protons to many com- 
mon impurities and the resulting ions may have large recombination coefficients. In this 
work, the mass spectrometer indicated that H 3 0 + ions formed by the fast (k= 4.3 x 10* 9 
cm 3 /s) proton transfer 


H 3 + + H 2 0 H 3 0 + + H 2 


( 8 ) 



were the principal offending ions. Hence, all efforts were made to reduce the H2O con- 
centration by placing cold traps into the gas supply lines, but some water vapor remained. 
Since the mass spectrometer is located far downstream from the recombination region, 
the mass spectra reflect conditions in the late afterglow, rather than in the recombination 
zone and some computations are required to find out how many H30 + ions are present in 
the recombination region. Model calculations were made which included proton transfer, 
recombination of both H3 + (using a= 2 xlO* 7 cm 3 /s ) and H30 + (a~lxl0‘ 6 cm^/s) [ 38 ] 
and ambipolar diffusion in the region from the H2 inlet to the mass spectrometer. The 
residual H2O concentration then was deduced by fitting the ratio H30 + /H3 + , as measured 
by the mass spectrometer, to that predicted by the model. This procedure yielded a 
residual water concentration of SxlO^cm' 3 . The model also shows that such a 
concentration of water will only have a small effect on the electron density decay in the 
recombination region (an additional 5% reduction during the first factor-of-ten decline of 
the electron density). 


II. 6 Data analysis 


Several different methods were used to derive de-ionization coefficients from the 
measured electron densities n e (z, r=0). The customary method of deducing P from the 

slope of a graph of 


1 

n e (z) 


1 

n e (z = 0) 



(9) 


was used to obtain approximate values of the de-ionization coefficient. In this simple 
approach, diffusion in axial and radial directions is neglected. A more rigorous approach 
is needed if one wishes to determine the dependence of P on z. It then becomes necessary 



to analyze the shape of the decay curves in detail and one is forced to include diffusion in 
both axial and radial directions. 

A rigorous analysis of a recombining flowing afterglow plasma is a rather difficult 
task. The continuity equation for the electron density n e (r, z) in the flow tube with a 
parabolic velocity field (Poiseuille flow) is given by 


dz 


1 




{ -k* d 4 


( 10 ) 


Here, V 0 is the gas flow velocity on the axis of the tube, and r and R denote the radial 
variable and the flow tube radius, respectively. D a is the ambipolar diffusion coefficient 
of ions and electrons. It is was deduced from measured mobilities using Einstein’s 
relation between mobility and diffusion. A series of measurements in H 3 + plasmas at very 
low electron densities (< 2xl0 9 cm' 3 ) showed that the calculated diffusion loss agreed 
with the measured diffusion loss to an accuracy of better than ± 10%. 

Numerical solutions of Eq. (10) were obtained for a few conditions of interest, 
but the computing time was found to be impractically long for fitting a large set of exper- 
imental data. Hence, the computed solutions were used only to validate a simpler, ap- 
proximate method that was employed for routine data analysis. In the approximate 
method, the gas flow velocity is taken to be radially independent and the assumption is 
made that the radial variation of n e (z, r) is given by fundamental diffusion mode in a 
cylinder. The approximate continuity equation for the electron density on axis (r=0) then 
can be written as 


dnj£ = 0 ) 
dz 


V 'ff 




(ID 



Here, A= R/ 2.405 is the fundamental diffusion length for a cylinder of diameter R and 
D a is the ambipolar diffusion coefficient.. The effective flow velocity V e ff is close to the 
gas velocity on the axis of the flow tube. Estimates showed that axial diffusion (the 
second-derivative term in Eq. 1 1 ) can have a significant effect in the early afterglow and 
that it should be retained even at this level of approximation. 

The z-dependent de-ionization coefficient P(z) was obtained by putting Eq. 1 1 
in the form 


1 


,d 2 n 




dn. 


( 12 ) 


If n e (z) is known from the measurement, all quantities on the right-hand-side are known 
and, at least in principle, p(z) can be computed. If only a single ion with recombination 
coefficient a is present, the procedure should return a constant value of P(z)=a. In 
practice, differentiating experimental data is a poor procedure and it was found 
preferable to fit the reciprocal values of n e (z) to a polynomial. The needed derivatives 
then were obtained analytically. A third-order polynomial fit to l/n e (z) results in a 
second-order polynomial expansion of P(z). Higher-order polynomials were tried, but 
they tended to introduce undesirable oscillations that were not present in the data. 

The approximate model was tested by using the exact model to generate simulated 
n e (z) values for known recombination coefficients and typical experimental conditions. 
The computer-generated data were then analyzed using Eq. 12. We found that the 
approximate method reproduced the input recombination coefficients (see Fig. 2 )to an 
accuracy of about 3% when the effective flow velocity was taken as 100% of the gas 
flow velocity on the axis of the flow tube. When the same method was applied to reduce 
experimental data obtained in recombining C>2 + plasmas (made by adding O 2 at the 
reagent inlet), the (3(z) curves exhibited only random deviations of <10% from an the 



average value of 1.85 x 10' 7 cm 3 /s. The average value of P(C> 2 + ) is only 6% smaller 
than the recombination coefficient given Mehr and Biondi [39] (a= 1.95 x 10' 7 cm 3 /s). 
This discrepancy may reflect, in part, a systematic error in the absolute values of n e (z) 
that are obtained from the Langmuir probes. 

During the data fitting it was noted that the P(z) curves at large z were very sensi- 
tive to small experimental errors in the n e (z) data at large z. This is not a flaw in the 
method of analysis. It merely reflects the fact that recombination becomes less important 
in the late afterglow. 

M. Experimental results 

A total of 40 decay curves n e (z) were recorded for different argon densities (from 
4x 10 15 cm' 3 to 8x 10 15 cm* 3 ) and different H 2 densities (from lx 10 14 cm -3 to 1.7x 
1015 cm* 3 ) . The absence of significant impurities was checked in each case by the mass 
spectrometer before the data were taken. 

The simplest form of data analysis consisted of deducing average de-ionization 
coefficients from the slope of l/n e (z) vs. z graphs. The results indicated that P increased 
noticeably with increasing H 2 concentration as is shown in Fig. 3. A similar graph was 
obtained when the average values of P(z) in the first 4 cm were plotted as a function of 
[H 2 ]. No such variation was found when the argon concentration was varied. 

Figures 4 and 5 show typical examples of the electron density n e (z) for two 
different H 2 concentrations. In the lower part of the figures, we show P(z) as calculated 
from Eq. 12. It is apparent that p decreases with increasing z. This finding qualitatively 
reproduces the observations of Smith and Spanel[13]. Actually, if their published data 
(Fig. 2 in ref.[13] ) are processed by using Eq. 12, the resulting P(z) curves look 
remarkably similar to those in our Fig. 4 and those shown by Smith and Spanel (their 
Fig.4 ). The decay of the electron density in this plasma clearly is anomalous in the sense 
that it cannot be described by a simple recombination law. 



Two comments should be made at this point: The (3(z) curves shown in Figs. 4 
and 5 may give the impression of higher accuracy than is justified by the original data. 
The initial values of (3(z) are reasonably insensitive to the choice of diffusion constants 
that enter in the analysis; the values at z=10 cm would be in error by ± 20 % if the 
diffusion losses have an error of ± 10 %. We also concede that the late-afterglow decay 
curves observed by Smith and Spanel [13] are intrinsically more accurate than ours. The 
reason is simply that the flow tube used by those authors had a larger cross section and 
that their experiment was done at a higher pressure. As a consequence, the diffusion 
losses were smaller (by a factor of 4) and the late afterglow data were superior. 

A smaller set of data was taken for D 3 + ions. The observations largely mirrored 
those for H 3 + . Again, the recombination coefficient increased with increasing D 2 concen- 
tration (see Fig. 3), but the D 3 + de-ionization coefficients were significantly smaller (from 
0.7xl0' 7 cm 3 /s to 1.3 xlO' 7 cm 3 /s in the range [D 2 ] = 1 to 14 x 10 14 cm' 3 ). 

IV. Search for excited states 

Several auxiliary measurements were carried out to find evidence for the 
presence of excited H 3 + ions. The ideas for those were, in part, based on the tests 
performed by Smith and Spanel. 

In one such test, carbon monoxide was added through a reagent inlet located 
about 6 cm downstream from the H 2 inlet. At that location, most of the fast recombining 
fraction of H 3 + should have recombined, but the slowly recombining fraction should still 
be present. The reaction with CO, 

H 3 + + CO — »HCO + + H 2 (13) 

would convert the non-recombining H 3 + ions to recombining HCO + (a = 2.4x 10' 7 
cm 3 /s) [40] ions. As a consequence, the electron density should fall of more rapidly after 



CO addition. Smith and Spanel observed a distinct enhancement of the electron loss upon 
adding CO. We tried to duplicate this experiment, but we found no effect; substitution of 
HCO + ions for H3 + ions did not measurably enhance the electron loss. We do not con- 
sider our failure to observe the effect of CO additions as being in conflict with Smith and 
Spanel's findings. The diffusion losses in our tube were higher, which makes it more 
difficult to observe subtle effects in the late afterglow. Addition of methane at the 
downstream location, however, did result in a clear enhancement of the electron loss 
rate. This is not surprising since the proton transfer with methane produces CHs + ions (a 
= l.lx 10 ' 7 cm 3 /s) [ 41 ] which recombine much faster than H3 + . The test with methane 
also shows that H3 + ions were still present at that point in the afterglow. If the plasma at 
that point had contained mainly H30 + impurity ions, adding methane would not have had 
the observed effect. 

An attempt was also made to use methane addition as a detector for H3 + (v>3), 
even though there are strong reason to doubt their presence. The idea again was borrowed 
from Smith and Spanel, who give more details. It relies on the fact that H3 + (v>3) reacts 
with CH4 to form CH3 + which then reacts again with CH4 to form C2H5 + . In contrast, 
H3 + (v<3) reacts with CH4 to form CH5 + . If CHs* and C2H5 + had different electron- 
ion recombination coefficients, this could be exploited as a means to diagnose the vibra- 
tional state distribution of the H3 + ions. As a preliminary step, we measured the recom- 
bination coefficients of both ions, however using Ar* ions reacting with CH4 to make 
CH3 + which then reacted further to make C2Hs + . To our dismay, we found only a 10 % 
difference between the recombination rates of CH5 + (a ~lxl0' 6 cmVs) and C2H5 + , 
while Smith and Spanel had found that CHs + recombined significantly faster. In view of 
this problem, we discontinued the search for H3 + (v>3) by this method. 

Smith and Spanel also described a rather interesting method of detecting H3 + ions 
in v=l which seemed worth repeating. They noted the presence of a small amount of 
ArH + ions in the downstream plasma. This was found to be surprising since all ArH + ions 



should have been destroyed in reactions with the abundant H 2 . Their explanation for the 
persistence of ArH + was that those ions entered an equilibrium with vibrationally excited 
H 3 + (v=1) ions via the reaction 

H 3 + (v=1) + Are ArH+ + H 2 . (14) 

The observation of ArH + then was seen as evidence for the existence of vibrationally ex- 
cited ions. We repeated this measurement in the hope that we might leam more about the 
abundance of H 3 + in v=l. To avoid complications due to recombination reactions of 
ArH + and H 3 + , we used our selected-ion drift tube apparatus [42 ] rather than the plasma 
flow tube. In this experiment, short pulses (10 |is) of mass selected Ar* ions were in- 
jected into an ion drift tube of length 36 cm, filled with helium buffer gas with additions 
of argon and H 2 . The reactions taking place are the same as those in the flow tube . A 
second mass spectrometer, located at the opposite end of the tube, analyzes the ions that 
have traversed the drift tube. The arrival time spectra of the ions were recorded with a 
time resolution of about 0.25 p.s. If two ion species, e.g. H 3 + (v=l) and ArH + , react 
back and forth many times while traversing the drift tube, their arrival spectra should be 
nearly identical. This was indeed observed. The ArH + signal, while far smaller (<1% of 
H 3 + ), exhibited to an extraordinary degree of precision the same time dependence as that 
of H 3 + . At first, this seemed to indicate an equilibrium process. However, it was 
discovered rather quickly that the ratio of the ArH + and H 3 + intensities did not vary 
linearly with the [Ar]/[H 2 ] concentration ratio, as would be required for an equilibrium 
process. Further tests showed that ArH + signal depended critically on the small bias 
voltages that were applied to the ion extraction electrodes and it became clear that the 
vast majority of ArH + ions were produced by proton transfer from non-thermal H 3 + ions 
during the extraction of ions from the drift tube into the mass spectrometer. We suspect, 
but cannot prove, that the same complication was present in the flow tube experiments of 



Smith and Spanel. The equilibrium "constants" given by them (actually they varied by a 
factor of four) may have little to do with the reaction equilibrium of reaction 14 and no 
conclusions about the existence of excited H 3 + ions can be drawn. 


V. Interpretation of Results 

The interpretation of the data should provide an explanation for the enhancement 
of the de-ionization coefficient due to hydrogen addition and the decline of the rate of 
recombination with position z. In the following section, we will examine two different 
models. The first model assumes the presence of H 3 + ions in two different states. The 
second model will invoke an entirely different explanation in terms of Rydberg 
molecules. 

V.l. Two-state model 

In view of our discussion in Sec. II, it is natural to seek an explanation of the z 
dependence of P by assuming that H 3 + ions in v= 0 and in vi=l are present. If this were 
the case, then the de-ionization coefficient for a mixture of v=0 and vj=l ions would be 
given by 


p (z) = oto f 0 (z) + ai fi(z) (15) 

where ao , oq denote the recombination coefficients of v=0 and vj=l ions, and fo and fi 
are the fractional concentrations (fo+ fi=l). Here the oq are constants, but the fractions 
f G and fi will change with z as a result of recombination and quenching of excited ions. 

A reaction model that reproduces the measured P(z) curve obtained at [H 2 ]= 1.3 

x 10 14 cm' 3 can be constructed easily. It consists of recombination of ions in v=0, vj=l, 
and quenching of vi=l ions to v=0 in collisions with H 2 . Four quantities enter: a G and 



dj, the initial fraction f o (z=0), and the quenching coefficient kq of the vi=l state with 
H 2 . An adequate fit (see Fig. 4 ) can be obtained by using a G = 2.2 xlO' 7 cm^s for the 
recombination coefficient of v=0 ions, aj= 4 xl0‘8 cm 3 /s for that of v=l ions, and f 0 = 
0.7. This choice of f Q was suggested by the results of Lee et al.[32]. Two trial values 
for the quenching coefficient were used, kq= lxlO * 17 cm 3 /s and a ten times larger value 
(it does not enter critically at low H 2 concentrations). As may be seen in Fig.4, the two- 
state model reproduces the measured de-ionization coefficients fairly and we initially 
believed that this model had some merit. 

A problem arises, however, when we use the same values to compute a fit to the 
data at a higher H 2 concentration, [H 2 ]= 1.7 x 10 15 cm -3 . The model calculations (see 
Fig. 5 ) indicate that p(z) should increase with z and remain nearly constant thereafter 
when the higher value of kq is used. This happens because quenching of the v^=l state 
now occurs much earlier in the afterglow. We can force a decline of p(z) by reducing 
the quenching rate to kq=lxlO' 13 cm 3 /s . If such a small value of kq is taken, however, it 
is hard to explain why the measured P at z=0 (only 1 cm from the reagent inlet) should be 
larger than at the lower H 2 concentration. We eventually concluded that the two-state 
model does not really provide a consistent explanation of the experimental findings. In 
the absence of independent evidence for the presence of excited states of H 3 + , it would 
not seem justified to introduce additional H 3 + states or other quenching mechanisms into 
the model in an attempt to explain the data. 

It is noteworthy that our value of p(z=0) = 2.2 xlO' 7 cm 3 /s at [H 2 ] = 1.7x 10 15 
cm 3 /s is very close to that by Leu et al. .[7] (2.3 xlO' 7 cm 3 /s). We should not have 
observed such high values if the plasma had contained a sizable fraction of non- 
recombining ions. The experiment of Leu et al. made use the very accurate microwave 
afterglow method which is free of many complications that are present in flowing 
plasmas. In retrospect. Leu et al.s' claim that their H 3 + ions were v=0 was well justified 



since its is now known that discharges in He/H 2 mixtures contain very few v=l ions 
(about 6% of v=0)[34], 

V.2 De-ionization by Rydberg molecule stabilization (RMS). 

The interpretation of our data in terms of two states of the H 3 + ion did not result 
in a consistent picture. If it had done so, we still would have had the problem of 
explaining how H 3 + ions in v=0 undergo recombination in the absence of suitable 
crossings of potential surfaces while those in vi=l do not. The task now will be to find a 
de-ionization mechanism that is compatible with the anomalous decay without invoking 
binary recombination. 

Our data contain a clue about the de-ionization mechanism. The slow decay of the 
plasma at late times (low electron densities) and the effect of H 2 suggest that both 
electrons and H 2 molecules in some way influence the process of recombination as such, 
rather than affecting the identity of the recombining ion species. Such considerations and 
information obtained from other experiments motivated us to explore an entirely different 
de-ionization mechanism. The process discussed in the following departs significantly 
from previous models of de-ionization in plasmas. In judging its merits, it should be kept 
in mind that no previous theory has provided a quantitative explanation of H 3 + 
recombination. 

The basic idea is that de-ionization of H 3 + plasmas occurs by a very efficient form 
of three-body recombination which involves both electrons and neutral atoms as third 
bodies. The dominant part of the recombination energy is disposed of in the form of 
radiation or by a reaction with H 2 . 

Three-body recombination, in which electrons act as stabilizing third bodies, is 
usually considered only in plasmas containing atomic ions (for a recent review of three- 
body recombination, see Flannery [43]). In this mechanism, the captured electron 
transfers part of its kinetic energy to a second electron while being accelerated in the ion's 



field. The process is quite well understood. It’s recombination coefficient [44 ] is given 
by 


a [cm 3 /s]= n e 2.5 x 1CH 9 (300/T) 9 /2 (16 ) 

which implies small values of a = 2.5 xlO’ 9 cm 3 /s at n e ~ 10 10 cm 3 /s and T=300 K. 
When molecular ions are present, dissociative recombination tends to be far more rapid 
and for this reason little attention has been paid to electron-stabilized recombination in 
plasmas that contain molecular ions. The same is true for neutral-stabilized 
recombinauon, which is less efficient by 6 to 7 orders of magnitude than electron- 
stabilized recombination. A special case arises, however, when a plasma contains 
molecular ions that do not undergo dissociadve recombination. Long-lived autoionizing 
states can then be formed in collisions of electrons with molecular ions, and 
redistribution of angular momenta (1-mixing) can further enhance the lifetimes of such 
states. 

The following sequence may occur in H3 + plasmas: In the first step, an electron 
of low angular momentum (s or possibly p) is captured into a v=l Rydberg state of H3 , 
i.e.. 


e* + H3+ — > H3** (v=l, n, s ). (17) 

There is nothing new in invoking this resonance. It is the first step in the indirect 
recombination mechanism [45]. Excitation to the V2=l state (0.3126 eV) in n= 7, 1=0, 
would be resonant at an electron energy of +0.0309 eV (1.2 kT at 300 K), very close to 
the maximum of a thermal distribution. For D3 + , capture into n=8 , 1=1 would be 
resonant at ~0.5 kT, which is somewhat less efficient. 



When there are no crossings between Rydberg states and repulsive states, 
predissociation will be rather slow and most doubly-excited states will decay by 
autoionization. If there were such crossings, then of course we would have a ready 
explanation for recombination, but this is not the case [3,5], 

It is now known from studies of dielectronic recombination [46] and zero- 
electron-kinetic-energy (ZEKE) spectroscopy [47 ] that redistribution of angular 
momenta of Rydberg electrons (1-mixing) due to electric microfields and collisions can 
reduce the rates of predissociation and autoionization by several orders of magnitude, the 
basic reason being that high-1 electrons do not penetrate to the core of the molecule. We 
propose that ambient plasma electrons, neutral atoms (He), and perhaps ions rapidly 
induce 1-mixing in the H3** molecules and render them stable (or long-lived) against 
autoionization. Since the statistical weight of states increases with 1, most H3** will not 
autoionize so that their average lifetime is drastically increased. A mixing of nearly 
degenerate 1-states and their magnetic sublevels requires only momentum exchange but 
very little energy exchange with third particles and thus should occur very rapidly. The 
expected mixing efficiency will be discussed next. 

If 1-mixing is to be effective, it must occur during the lifetime of the autoionizing 
resonance. The lifetime of some H3** autoionizing states have been measured to be 
extremely long (of the order of 10*7 s)[48]. In particular, autoionizing states formed by 
e *+ H3 + collisions in an ostensibly field-free vacuum region must survive for 10*8 to 10*? 
s to explain observations in the merged-beam experiment of Yousif et al. [21] and those 
of Mitchell et al. [49] . Unfortunately, the state of the long-lived H3** molecules that are 
observed in the beam experiments, is not well known. It seems possible that the electron 
is captured into a state of low angular momentum and then exchanges angular momentum 
with the ion-core. Perhaps, this may be the reason that the lifetimes are so long. It is also 
not entirely clear whether or not the observed lifetimes really refer to molecules that have 
not interacted in any way with the small electric stray fields (estimated by Mitchell [49] 



to be < lV/cm) that are always present in beam experiments. Perhaps, some 1-mixing 
occurs even in small, static fields of <1 V/cm. If this were the case, however, then we 
would argue that 1-mixing in a plasma should be at least as efficient but probably far 
more efficient. The electric microfields in a plasma at n e =2xl0^ cm*3 due to nearest- 
neighbor charge (distance of a few pm) are also of the order of lV/cm and the fields have 
high-frequency components up to 10 11 Hz. 

Measurements of Gallagher et al.[50] on sodium Rydberg atoms in n= 5-10 
indicate that 1-mixing in collisions with helium atoms proceeds with a rate coefficient of 
~ 10' 8 cm 3 /s if 1>2. At 1 Torr of He, this gives a mixing time of 3 ns which is shorter than 
the lifetime of the H3** resonance. However, for 1=0 and 1, mixing due to heavy- 
particle perturbers is not observed. The reason for the lack of 1-mixing of the s and p 
states is that those states have larger quantum defects (lower energy), and that heavy 
particles can transfer only a minute fraction of their energy to electrons. The quantum 
defects for H3 Rydberg states have been calculated by King and Morokuma [51] to be 
approximately 0.06 (s-states) and 0.025 (p-states), which means that an n=7 s-state would 
lie 0.005 eV=kT/5 below states of high 1 values. The transfer from s to p and higher 
states is the point where 1-mixing due to electrons may play a crucial and rate- 
controlling role. The energy barrier of kT/5 is no hindrance for thermal electron collisions 
and a single electron collision can transfer sufficient energy. Estimates of 1-changing 
cross sections [52] for 25eV electrons with n=25 atoms indicate very large values of 
10' 10 to 10‘ 9 cm 2 . Of more direct relevance are the experiments by Devos et al. [53] . 
They showed that the 1-changing efficiency of thermal electrons in collisions with helium 
Rydberg atoms (n= 8 to 17) is so large that mixing of 1-states occurs in an immeasurable 
time (less than a few ns) at electron densities of 2 xlO 11 cm‘3 . This implies an 1-changing 
rate coefficient for thermal electrons of >10'^ cnvVs! An estimate of the three-body 
coefficient for forming H3** in high 1-states from H3 + and two electrons would be given 
by the product of the capture rate coefficient (~a few 10‘ 2 cmVs), the lifetime of the 



resonance (~10' 8 to 10' 7 s), and the 1-mixing coefficient due to electrons (>10' 3 cm 3 /s). 
The resulting rate coefficient (~10' 17 cm 6 /s) would be larger that in Eq. 16 by two orders 
of magnitude. 

We will now turn our attention to decay modes of H 3 **(n=7, v=l, 1 ) molecules 
other than autoionization. They can be re-ionized by slow electron collisions which 
change their 1-value back to the autoionizing state (s, or p). In addition , they may 
descend to lower energy states by collisional energy loss with neutral atoms or electrons, 
slow predissociation (if it occurs), by radiating to a lower state, or by reacting with an H2 
molecule. Collisional energy loss (n-changing) with helium atoms is not important at 
n=7 (Bates and Khare [54]) compared to radiative decay. The rate coefficient for n- 
changing collisions of electrons with n=7 Rydberg atoms (electron quenching) has been 
measured by Devos et al. [52 ] to be ~10' 5 cm 3 /s which yields a destruction rate of 
2x10^ s'l at n e = 2 xlO 1 * 3 cm' 3 . Radiative decay rates for n=7 hydro genic states are 
close to lxlO 6 s' 1 for s-states , 2x 10 7 s' 1 for p-states, and about 2xl0 6 s' 1 for a 
statistical mixture of 1- values up to 1=6 [51, 55 ] . Therefore, radiation is likely to be the 
dominant decay mode under the conditions of our experiment, but electron quenching 
will dominate at higher electron densities. Destruction in collisions with H2 will be a two 
step process, formation of an H5* complex, followed by predissociation, 

H 3 ** + H 2 H 5 * H 2 + H 3 (n=3, 4 ). (18) 

The second step is essentially the same as dissociative recombination of H5 + ions which 
is known to occur rapidly. Reaction (18) is important for the interpretation of our 
experiment, since it provides the needed enhancement of the loss rate of electrons in the 
presence of H2. The reaction will compete with radiative decay if it occurs on a time 
scale of clfis . At [H2] = lx lO 1 ^ cm' 3 , this requires a rate coefficient of ~ 1x10'^ cm 3 /s 
for reaction 1 8 , which is not unreasonable. 



In the following, we present a simple rate model for the de-ionization process 
and compare it to our data. The capture and autoionization step 

H 3 + + e - <-» H 3 ** (n=7, s,p) (19) 

is described by a capture coefficient kc and an autoionization rate v a . Stabilization by 1- 
mixing due to electrons 

H 3 ** (n=7, s,p)+ e ' <-> H 3 ** (n=7, l>l)+e * (20) 

and its reverse reaction are described by rate coefficients k s and kj in the forward and 
reverse directions, respectively. Helium does not enter as a third body, since we have 
assumed that the electrons control the rate-limiting step in the 1-mixing from s or p to 
higher 1 and back.The main destruction reactions for H 3 ** are: 


H 3 **(n=7) -> 

H 3 * (n<7) + hv 

(radiation, rate v,^ ) 

(21a) 

H 3 ** -> 

H 2 + H ( or 3H ) 

(predissociation, rate v p ) 

(21b) 

H 3 **(n=7) + e- -» 

H 3 *(n<7) +e* 

(rate coefficient kq e ) 

(21c) 

H 3 ** + H 2 — » H 2 + H 3 (low n) 

( rate coefficient kn ) . 

(2 Id) 


In order to allow for the possibility that the rate coefficients for H 3 ** in the 
autoionizing state (s,p) differ from those in higher 1-states, we treat their destruction 
separately and attach the subscripts "0" and "1" to the respective rate coefficients. It is a 
simple matter to show that this reaction scheme gives a de-ionization coefficient 
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where 

Vo = Vrad o + V p> o + kq e + k H [ H 2 ] 
and Vi = Vrad,! + v p> i + kq e tig + k H [ H 2 ]. 


( 22 ) 


(23a) 

(23b) 


In the limit of large n e and high values of v, (3 approaches the constant value k^ This 
means that all captured electrons recombine irreversibly. In the limit of vanishing n e , the 
deionization rate coefficient (3 approaches the value 

p — V- < 24 > 

1+2-s- 

v 0 

This value still depends on the H 2 concentration through Eq. 23a. In the limit of 
vanishing [H 2 ], (3 should approach a small value that depends on the ratio Va/v ra d,0- The 
ratio may be on the order of ten to hundred and (3 would then would be 1 to 10 % of kc. 

V.3 Comparison of the RMS model to experimental data 

The functional form of Eq. 22 fits our data quite well, and it is possible to find a 
set of rate coefficients such that Eq. 22 provides a reasonable fit. Fig. 6 shows a 
comparison to experimental values of p (n e ) (the same data as shown in Fig. 4 and 5, but 
plotted as a function of the electron density). The rate coefficients used were: 
kc = 4.4 x 10* 7 cm 3 /s 
v a = 1.3xl0 8 s' 1 



k s = 1.25 xlO' 2 cm-Vs 
ki =3.1 xl(H cm 3 /s 
Vrad,0= 1 - lx 10 6 s' 1 
Vrad,l= 1-1 x 1 0 6 s - 1 

Vp,0=Vp,l =0 

kn = 4.5 x 10‘ 9 cm^/s 

kq e = lx 10' 5 cm 3 /s 

This set of numbers is not necessarily unique but the values are quite close to 
those that were estimated earlier. Detailed balancing suggests that the ratio ks/kj should 
be close to n 2 -l=48 (assuming that only the 1=0 state autoionizes), where n=7 is the 
principal quantum number of the Rydberg state, since 1-mixing from s to 1>0 states 
requires little energy exchange(0.2 kT). The ratio ks/kj =40 fits the data quite well. 

Our model is crude in many respects. For instance, we treat all Rydberg 
molecules with 1>1 as a single entity which is probably not justified. We also do not wish 
to give the impression that our experimental data provide sufficient evidence for the 
validity of this model. More detailed studies would be needed to do that. On the other 
hand, the estimates underlying the model seem sound and the rate coefficients are 
reasonable. 

V. 4 Consequences of the proposed de-ionization mechanism and relation to other 
work. 

If the proposed de-ionization mechanism occurred, it would immediately solve 
the theoretical dilemma of having to account for "binary recombination" of electrons with 
H3 + ions. This process is not needed to explain the observed de-ionization. It is only 
required that the H3** Rydberg molecules decay to lower states by radiation, slow 



predissociation, or collisional quenching, and this will certainly occur. Radiation alone 
would suffice to explain the observed rate coefficients. 

In support of our model, we adduce several pieces of supporting experimental 
evidence. The first of those is somewhat indirect but it shows that "loosely bound" 
electrons (Rydberg molecules) can persist in some plasmas. The presence of "loosely 
bound electrons" in pure-helium afterglows has actually been a puzzle for many years. 
He2 + ions dominate such plasmas, and they are known noi to undergo rapid dissociative 
recombination. The situation may be quite similar to that in H3 + plasmas. Experiments by 
Kaplafka et al. [ 56 ] showed that exposure of a helium afterglow plasma to a rather weak 
pulse of microwave radiation (capable of heating free electrons by several hundred K) 
enhanced the free electron density by factors of up to three. This would indicate that the 
plasma contained more loosely-bound electrons than free electrons. Similar observations 
have been made by Castell [ 57 ] . It seems very likely that the "loosely bound electrons" 
are He2* Rydberg molecules formed by collisional stabilization of He2** autoionizing 
complexes, the same process that we invoke in the case of H3 + plasmas. It can be shown 
from the equations given earlier that the concentration of Rydberg molecules can 
approach or exceed that of free electrons if the electron density is sufficiently high. 

To our knowledge, no search for loosely bound electrons has been carried out in 
plasmas containing H3 + ions. There is, however, one puzzling spectroscopic observation 
that could be explained readily by postulating the presence of Rydberg molecules. Both 
Amano [ 58 , 59 ] and also Miderski and Gellene [ 60 ] noted that H3 (or D3) absorption and 
emission lines (from n =3 states) in H2 (D2) plasmas exhibited line profiles that indicated 
high, non-thermal energies (about 0.3 eV) of the absorbing or emitting H3 molecules. As 
noted by the authors, the energy is close to that expected if the excited (n= 3 ) H3 
molecules were formed by recombination of H5 + . But, as Amano noted, the plasma 
should not have contained any H5 + ions under his experimental conditions. The solution 
to the seeming contradiction may be that the fast, excited H3 species were produced from 
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H3* + H2 collisions, Eq. 18 , which has nearly the same energetics as H5 + 
recombination. Thus, spectroscopic observations may give the appearance of Hs* 
recombination even though these ions are not present in the plasma. Actually, there may 
be two excitation mechanisms for H3(n) states, decay from H3** (no excess translational 
energy) and reaction ( 18 ). At high lie and low T, collisional radiative recombination (Eq. 
16 ) should provide an additional mechanism of population Rydberg states. Spectroscopic 
studies may be able to distinguish the different excitation mechanisms. 

The most direct observation that indicates production of H3** molecules from 
H3 + recombination comes from the merged-beam experiments of Mitchell [ 49 l. He found 
that the apparent H3 + recombination cross sections could be increased by about a factor 
of five when a "deflection field" (used to deflect unrecombined ions out of the beam 
before they strike the detector) was reduced in magnitude. This would suggest that 
majority of H3 + + e' "recombination" events produces long-lived H3** Rydberg 
molecules. The role of electric fields in beam experiments may be quite complicated and 
probably deserves some more analysis. 

The proposed mechanism, if it is valid, has also one disturbing consequence: 
Observation of plasma de-ionization coefficients on the order of a few 10‘ 7 cmVs cannot 
be taken as sufficient evidence for a dissociative mechanism. Most ions that have been 
studied in the past probably do recombine dissociatively, but there is very little direct 
evidence in support of that assumption. Stabilization of Rydberg molecules is likely to be 
more important in polyatomic ions and in may, in part, explain the large recombination 
rates of those. It may also play a role as an additional recombination channel in some 
simpler ions, possibly HCO + , HeH + and other hydride ions, and it may also give rise to 
some reaction products that cannot be produced by binary recombination. Thus, it can 
probably not be taken for granted that the product yields measured in afterglow 
experiments are applicable to environmejjftnts at low densities, such as many 
astrophysical plasmas. 
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VI. Conclusions 

Our experimental observations, in part, are very similar to those of Smith and 
Spanel [13 ] and we concur with them that the decay of H3 + plasmas is indeed 
anomalous in the sense that it does not obey a simple recombination law. However, our 
interpretation of the experimental findings is entirely different. 

We attempted to explain the anomalous plasma decay by assuming the presence 
of ions in the v=0 and vj = l state, but we did not find that this resulted in a consistent 
description of the experimental data. An alternative de-ionization model is proposed that 
invokes collisional stabilization of H3** Rydberg molecules. If this model is correct, then 
it is quite possible and even likely that the binary recombination of electrons with H3 + 
ions is a very slow process, as was predicted by the early theory of Michels and Hobbs 
[5]. 

Our measurements indicate that de-ionization coefficients in the range from 1.5 
to 2.5 xlO' 7 cm 3 /s at 300 K may be appropriate for modeling of H3 + plasmas of 
reasonably high density, for instance those in electrical discharges and perhaps in some 
planetary ionospheres. The T -0 - 5 dependence of oc(H3 + ) on gas temperature that was 
measured by Leu et al.[7] may be of some practical value to scale this "recombination 
coefficient" to other gas temperatures. The de-ionization coefficient of D3 + appears to be 
significantly smaller than that of H3 + . 

It is an enurely different quesdon which recombinauon coefficient should be used 
in environments of very low densities, e.g. in the interstellar medium. Here, the "true" 
binary recombination coefficient at low temperatures is needed. It is conceivable that the 
recombination cross sections measured in ion-storage ring experiments are close to the 
true values, but this is far from obvious since the presence of electric fields in the 
interaction region may also lead to 1-mixing effects. The recent observations by Mitchell 
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and coworkers [49] indicate that the interpretation of beam experiments may be far more 
complicated than had been thought. 

Acknowledgments: 

This work was, in part, supported by NASA under Grant No. NAGW 1764 

The authors gratefully acknowledge the assistance of A. Prabhaharan in carrying out the 
drift tube experiments, and the help of E.V. Shun'ko in the design and testing of the 
Langmuir probes. One of us (R.J.) acknowledges several helpful discussions with J.B.A. 
Mitchell, N.G. Adams, R. Flannery, and the late Sir David Bates. 

References: 

1 A. Dalgamo, in Advances in Atomic, Molecular and Optical Physics, B. 

Bederson and H. Walthers eds., Vol. 32, 57 (1993) 

2 T. Oka and M.-F. Jagod, J. Chem. Soc. Faraday Trans. 89, 2147 (1993) 

3 D. R. Bates, M.F. Guest and R.A. Kendall, Planet. Space Sci. 41, 9 (1993) 

4 J.N. Bardsley, J. Phys. B 1, 365 (1968) 

5 H.H. Michels and R.H. Hobbs, Ap. J. (Letters) 286, L27 (1984) 

6 K.C. Kulander and M.F. Guest, J. Phys. B 12, L501 (1979) 

7 M.T. Leu, M.A. Biondi and R. Johnsen, Phys. Rev. A 8, 413(1973) 

8 J. A. Macdonald, M.A. Biondi and R. Johnsen, Planet. Space Sci. 32, 651 (1984) 

9 N.G. Adams, D. Smith and E. Alge, J. Chem. Phys. 81, 1778 (1984) 


10 


R. Johnsen, Int. J. Mass Spectrom. Ion Processes 81, 67 (1987) 



34 


11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 


23 


A. Canosa, J.C. Gomet, B.R. Rowe, J.B.A. Mitchell, and J.L. Queffelec, J. Chem. 
Phys. 97, 1028 (1992) 

A. Canosa, B.R. Rowe, J.B.A. Mitchell, J.C. Gomet, and C. Brion, Astron. 
Astrophys. 248, L 19 (1991) 

D. Smith and P. Spanel, Int. J. Mass Spectrometry and Ion Processes 129 163 
(1993) 

R. Johnsen, A. Chen, and M.A. Biondi, J. Chem. Phys. 72, 3085 (1980) 

T. Amano, Ap. J. 329, L121 (1988) 

T. Amano, J. Chem. Phys. 92, 6492 (1990) 

M. Feher, A. Rohrbacher, J.P. Maier, Chem. Phys. 185, 357 (1994) 

H. Hus, F.B. Yousif, A. Sen, and J.B. A. Mitchell, Phys. Rev. A 38, 658 (1988) 

J.B.A. Mitchell, J.L. Forland, C.T. Ng, D.P. Levac, R.E. Mitchell, P.M. Mul, W. 
Claeys, A. Sen, and J. Wm. McGowan, Phys. Rev. Lett. 51, 885 (1983) 

J.B.A. Mitchell, C.T. Ng, J.L. Forland, .R. Janssen, and J. Wm. McGowan, J. 
Phys. B 17, L909 (1984) 

F. B. Yousif, P.J.T. Van der Donk, M. Orakzai, and J.B. A. Mitchell, Phys. Rev. A 
44 , 5653 (1991) 

M. Larsson, H. Danared, J.R. Mowat, P. Sigray, G. Sundstrom, L. Brostrom, A. 
Filevich, A. Kallberg, S. Mannervik, K.G. Rensfelt, and S. Datz, Phys. Rev. Lett. 
70 , 430(1993) 

G. Sundstrom, J.R. Mowat, H. Danared, S. Datz, L. Brostrom, A. Filevich, A. 
Kallberg, S. Mannervik, K.G. Rensfelt, P. Sigray, M. af Ugglas, and M. Larsson, 
Science 263,785,(1994) 



35 


24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


B. M. Dinelli, S. Miller, and J. Tennyson, J. Mol. Spectr. 153, 718 (1992) 

S. L. Guberman, Phys. Rev. A 49, R4277 (1994) 

D. R. Bates, Mon. Not. R. Astron. Soc. 263 , 369 (1993) and D. R. Bates, Proc. 
Roy. Soc. Lond. A 443 , 257 (1993) 

R. Johnsen, E.V. Shun'ko, T. Gougousi, and M.F. Golde, Phys. Rev E (in print) 

A.A. Radzig and B.M. Smirnov, Reference data on atoms, molecules, and ions. 
(Springer Verlag 1980) 

D.L. Albritton, Atomic Data and Nuclear Tables 22, 1 (1978) 

G. C. Lie and D. Frye, J. Chem. Phys. 96, 6784 (1992) 

T. Amano, in "Dissociative Recom bination :Theory, Experiment and 
Applications", edited by B.R, Rowe, J.B. Mitchell, and A. Canosa (Plenum Press, 
New York 1993) 

H. S. Lee, M. Drucker, N.G. Adams, Int. J. Mass Spectrom. Ion Processes 1 17, 
101 (1992) 

C. R. Blakley, M.L. Vestal, and J.H. Futrell, J. Chem. Phys. 66, 2392 (1977) 

M.G. Bawendi, B.D. Rehfuss, and T. Oka, J. Chem. Phys. 93, 6200 (1990) 

S. L. Guberman, in "Dissociative Recombination:Theory, Experiment and 
Applications", edited by B.R, Rowe, J.B. Mitchell, and A. Canosa (Plenum Press, 
New York 1993) 

T. Oka, Philos. Trans. R. Soc. London, A, 303 , 543 (1981) 

R. Johnsen, C.-M. Huang, and M.A. Biondi, J . Chem. Phys. 65, 1539 (1976) 

M.T. Leu, M.A. Biondi and R. Johnsen, Phys. Rev. A 7, 292 (1973) 

F.J. Mehr and M.A. Biondi, Phys. Rev. 181, 264 (1969) 



40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


B. Ganguli, M.A. Biondi, and R. Johnsen, Phys. Rev. A 37, 2543 (1988) 

N.G. Adams, D. Smith, and E. Alge, J. Chem. Phys. 81, 1778 (1984) 

R.E. Tosh and R. Johnsen, Int. J. Mass Spectr. and Ion Proc. 123, 193 (1993) 

M.R. Flannery, Advances in Atomic, Molecular and Optical Physics 32, 1 17 
(1994) 

H.S.W. Massey and H.B. Gilbody, Electronis and Ionic Impact Phenomena, Vol 
4, page 2166 ( Oxford at the Clarendon Press, 1974) 

J.N. Bardsley and M.A. Biondi, Advances in Atomic, Molecular and Optical 
Physics, Vol. 6, 1 (1970) 

A. Muller, D.S. Belie', B.D. DePaola, N. Djuric', G.H. Dunn, D.W. Mueller, and 

C. Timmer, Phys. Rev. A 36, 599 (1987) 

F. Merkt and R.N. Zare, J. Chem. Phys. 101, 3495 (1994) 

V. Berardi, N. Spinelli, R. Velotta, M. Armenante, and A. Zecca, Phys. Rev. A, 
47, 986 ((1993) 

J.B.A. Mitchell, Bulletin of the American Physical Society Series II, 39 1456 
(1994), and private communication 

T.F. Gallagher, S.A. Edelstein, and R.M. Hill, Phys. Rev. Lett. 35, 644 (1975) 

H.F. King and K. Morokuma, J. Chem. Phys. 71, 3216 (1979) 

T. F. Gallagher, "Rydberg Atoms", page 285, Cambridge University Press 1994 

F. Devos, J. Boulmer, J.-F. Delpech, J. de Physique 40, 215 (1979), 
and Phys. Rev. Lett. 39, 1400 (1977) 

D. R. Bates and S.P. Khare, Proc. Phys. Soc. 85, 231 (1965) 

J. R. Hiskes and C.B. Tarter (Report UCRL-7088, Lawrence Radiation 
Laboratory, Livermore CA, (1964)) 

J.P. Kaplafka, H. Merkelo, and L. Goldstein, Phys. Rev. Lett. 14, 970 (1968) 

F.R. Castell, Doctoral Thesis, University of Pittsburgh (1973) 


T. Amano, in "Dissociative Recombination:Theory, Experiment and 
Applications", edited by B.R, Rowe, J.B. Mitchell, and A. Canosa (Plenum Press, 
New York 1993) 


59 


T. Amano, SPIE Vol. 1858 (1993) 



60 


C.A. Miderski and G.I. Gellene, J. Chem. Phys. 88, 5331 (1988) 



Table I 

Recombination coefficients observed in earlier experiments. The abbreviations of 
methods refer to: SA (Stationary Afterglow), FALP (Flowing Afterglow/Langmuir 
Probe), IR (Infrared absorption in afterglow plasmas), MB (Merged Beam), ISR (Ion 
Storage Ring). The cross sections measured in beam experiments (MB and ISR) have 
been converted to 300 K rate coefficients. 

T Method vibrational reference 

[IQ- 7 cm 3 /s] [K] _state 

2,3±0,3 300 _SA 0 7 

1.6 _300 _SA _0 _8 

<0.2 _300 FALP 0 9 

1.8±0.2 210 IR _0 15 

1.5 (1.1) 300(650) FALP <3 (0) 11, 12 

<2 _510 JR _0 17 

1.3 (0.3) 300 FALP >3 (1,2) 13 

02 (300) MB 0 _21 

1.15 (300) ISR 0 22 
















































Figure Captions 

Fig-1 Schematic diagram of the flow tube used in these experiments. The main section 
of the tube has a length of 34 cm and an internal diameter of 6 cm. 

Fig. 2 Values of P(z) (filled squares, right-hand scale) derived from simulated ng(z) data 

(open circles, left-hand scale) by the procedure described in the text. Line drawn through 
the n e (z) data represents a third order polynomial fit. A recombination coefficient of 
2xl0' 7 cm 37 s was used in the simulation. The arrow indicates the position where the 
reagent is added in the simulation. 

Fig- 3 Measured de-ionization coefficients P in H 3 + and D 3 + plasmas obtained from the 
l/n e vs. z plots, as a function of H 2 (D 2 ) concentration. Heavy line: Least-squares fit to 
the data. Dotted lines represent the 95% confidence limits of the least-square fit. 

Fig. 4 Values of P(z) (filled squares, right-hand scale) derived from measured n e (z) data 
(open circles, left-hand scale) at [H 2 ] = 1.3xl0 14 cm -3 . Line drawn through the ne(z) data 
represents a third-order polynomial fit. The arrow indicates the position where the 
reagent is added in the simulation. Lines through P(z) points represent model 
calculations using the two-state model for a quenching coefficient of kq= lx 10' 12 cm 3 /s 
(dashed line) and kq= lxlO- 13 cm 3 /s (dotted line). 

Fig. 5 Values of P(z) (filled squares, right-hand scale) derived from measured n e (z) data 
(open circles, left-hand scale) at [H 2 ] = 1.7xl0 15 cm' 3 . Line drawn through the n e (z) data 
represents a third-order polynomial fit. The arrow indicates the position where the 
reagent is added in the simulation. Lines through P(z) points represent model 
calculations using the two-state model for a quenching coefficient of kq= lx 10' 32 cm 3 /s 
(dashed line) and kq= lxlO- 13 cm 3 /s (dotted line). 

Fig. 6 Comparison of observed values of P(n e ) at [H 2 ] = 1.3xl0 14 cm* 3 (circles) and 
[H 2 ] = 1.7xl0 15 cm* 3 ( squares) to the model calculations described in the text. 
















